Electron screening and its effects on Big-Bang nucleosynthesis 
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We study the effects of electron screening on nuclear reaction rates occurring during the Big Bang 
nucleosynthesis epoch. The sensitivity of the predicted elemental abundances on electron screening 
is studied in details. It is shown that electron screening does not produce noticeable results in the 
abundances unless the traditional Debye-Hiickel model for the treatment of electron screening in 
stellar environments is enhanced by several orders of magnitude. The present work rules out electron 
screening as a relevant ingredient to Big Bang nucleosynthesis confirming a previous study by Itoh 
et al. [l| and ruling out exotic possibilities for the treatment of screening, beyond the mean-field 
theoretical approach. 
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During the Big Bang the Universe evolved very rapidly 
and only the lightest nuclides (e.g., D, 3 He, 4 He, and 
7 Li) could be synthesized. The abundances of these nu- 
clides are probes of the conditions of the Universe dur- 
ing the very early stages of its evolution, i.e., the first 
few minutes. The conditions during the Big Bang nu- 
cleosynthesis (BBN) are believed to be well described 
in terms of standard models of cosmology and particle 
physics which determine the values of, e.g., tempera- 
ture, nucleon density, expansion rate, neutrino content, 
neutrino-antineutrino asymmetry, etc. Deviations from 
the BBN test the parameters of these models, and con- 
strain nonstandard physics or cosmology that may al- 
ter the conditions during BBN 0, H[ . Sensitivity to the 
several parameters and physics input in the BBN model 
have been investigated thoroughly in the past (see, e.g., 
) ■ In this work we consider if the screening by elec- 
trons would have any impact on the BBN predictions. 
This work reinforces the conclusions presented in Ref. 
PH , namely, that screening is not a relevant ingredient of 
the BBN nucleosynthesis. 

Modeling the BBN and stellar evolution requires that 
one includes the information on nuclear reaction rates 
(av) in reaction network calculations, where a is the 
nuclear fusion cross sections and v is the relative ve- 
locity between the participant nuclides. Whereas v is 
well described by a Maxwell-Boltzmann velocity distri- 
bution for a given temperature T, the cross section a 
is taken from laboratory experiments on earth, some of 
which are not as well known as desired ka-|8j]. The pres- 
ence of atomic electrons in laboratory experiments also 
influence the measured values of the cross sections in a 
rather unexpected way (see, e.g. [l(| )■ In the network 
calculations for the description of elemental synthesis in 
the BBN or in stellar evolution one needs account for the 
differences of the "bare" nuclear cross sections obtained 
in laboratory measurements, ab(E), to the corresponding 
quantities in stellar interiors, a s (E). One of these correc- 
tions is due to stellar electron screening as light nuclei in 
the stellar environments are almost completely ionized. 



Using the Debye-Hiickel model, Salpeter showed 
that stellar electron screening enhances cross sections, 
reducing the Coulomb barrier that reacting ions must 
overcome, yielding an enhancement factor 



f(E) = 



<t s (E) 
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The Debye-Hueckel model used by Salpeter yields a 
screened Coulomb potential, valid when (V) <C kT (weak 
screening) , 



(2) 



which depends on the ratio of the Coulomb potential at 
the Debye radius Rd to the temperature, 
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is the Debye radius, n is the ion number density, p is 
dimensionless quantity measured in units of g/cm 3 , 



7 2 

i>4 



1/2 



(5) 



where Xi is the mass fraction of nuclei of type i, and Tq 
is the dimensionless temperature in units of 10 6 K. The 
factor x corrects for the effects of electron degeneracy 

n 

Corrections to the Salpeter formula are expected at 
some level. Nonadiabatic effects have been suggested as 
one source, e.g., a high Gamow energy allows reacting 
nuclei having velocities significantly higher than the typ- 
ical ion velocity, so that the response of slower plasma 
ions might be suppressed. Dynamic corrections were first 

discussed by [l2| and later studied by Subsequent 

work showed that Salpeter's formula would be valid inde- 
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thermodynamic equilibrium of the solar plasma [lil - flfjj ]. 
Later, a number of contradictions were pointed out in in- 
vestigations claiming larger corrections, and a field the- 
oretic aproach was shown to lead to the expectation of 
only small (~ 4%) corrections to the standard formula, 
for solar conditions fl7l ]. 

Controversies about the magnitude of the screening ef- 
fect have not entirely died out and continued in some 
works [l8l - [2l| . These works are invariably based on 
molecular dynamics simulations. Dynamic screening be- 
comes important because the nuclei in a plasma are much 
slower than the electrons and are not able to rearrange 
themselves as quickly around faster moving ions. Since 
nuclear reactions require energies several times the av- 
erage thermal energy, the ions that are able to engage 
in nuclear reactions in the stars are such faster mov- 
ing ions, which therefore may not be accompanied by 
their full screening cloud. In fact, dynamic effects are 
important when particles react with large relative ve- 
locities. It appears that pairs of ions with greater rel- 
ative velocities experience less screening than pairs of 



particles with lower relative velocities [18j|. The correc- 
tion for dynamical screening can be approximated by re- 
placing Rd in eq. ([3]) by a velocity dependent quantity 

Rd(v p ) — Rd\Ji + fiVp/kT, where fj, is the ion pair re- 
duced mass and v p their relative velocity [22j . This effect 
is relevant in stellar environments whenever nuclear reac- 
tions occur at energies that are greater than the thermal 
energy. 

Experimentalists have exploited surrogate environ- 
ments to test our understanding of plasma screening ef- 
fects. For example, screening in d(d,p)t has been stud- 
ied for gaseous targets and for deuterated metals, insula- 
tors, and semiconductors (23l . |24| . It is believed that the 
quasi-free valence electrons in metals create a screening 
environment quite similar to that found in stellar plas- 
mas. The experiments in metals seem to have confirmed 
important predictions of the Debye model, such as the 
temperature dependence U e (T) oc T -1 / 2 24]. But there 
are still controversies on the validity of the experimental 
analysis and the use of the Debye screening, or Salpeter 
formula, to describe the experiments [III [26[ . 

A good measure of the screening effect is given by the 
screening parameter given by T — ZiZ%e 2 / (r)kT , where 
(r) = n -1 / 3 . In the core of the sun densities are of the or- 
der of p ~ 150 g/cm 3 with temperatures of T ~ 1.5 x 10 7 
K. For pp reactions in the sun, we thus get T ~ 1.06 
which validates the weak screening approximation, but 
for p 7 Be reactions one gets r p 7 Bo ~ 1.5, which is one of 
the reasons to support modifications of the Salpeter for- 
mula. Also, in the sun the number of ions within a sphere 
of radius Rd (Debye sphere) is of the order of N ~ 4. As 
the Debye-Hueckel approximation is based on the mean 
field approximation, i.e., for N = n(4irR 3 /3) 3> 1, devi- 
ations from the Salpeter approximation are justifiable. 

Based on the discussion above, there is a possibility 
that the screening enhancement factor, eq. ((T|), could 
appreciably differ from the Salpeter formula under sev- 
eral circumstances, leading to non-negligible changes in 
the BBN and stellar evolution predictions. It is the pur- 




T (10 9 K) 

FIG. 1: (Color online) Number density of electrons (solid 
curve) and of baryons (dashed curve) during the early universe 
as a function of the temperature in units of billion degrees 
Kelvin, Tg. 



pose of this work to verify under what conditions would 
this statement be true. 

In this work, the BBN abundances were calculated 
with a modified version of the standard BBN code de- 
rived from Wagoner, Fowler, and Hoyle [3j and Kawano 
[13, HI] (for a public code, see (29j). 

The electron density during the early universe varies 
strongly with the temperature as seen in figure [TJ where 
Tg is the temperature in units of 10 9 K (Tg). This can 
compared with the electron number density at the cen- 
ter of the sun, n™ ~ 10 26 /cm 3 . The figure shows that, 
at typical temperatures Tg ~ 0.1 — 1 during the BBN 
the universe had electron densities which are much larger 
that the electron density in the sun. However, in contrast 
to the sun, the baryon density in the early universe is 
much smaller than the electron density. The large elec- 
tron density is due to the e + e _ production by the abun- 
dant photons during the BBN. 

The baryonic density is best seen in figure [2] It varies 

as 



Pb 
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where h is the baryon density parameter 13011 . It can 
be calculated by using Eq. (3.11) of Ref. [3CJ and the 
baryon-to-photon ratio rj = 6.19 x 10~ 10 at the BBN 
epoch (from WMAP data (HI). Around Tg ~ 2 there 
is a change of the value of h from h ~ 2.1 x 10~ 5 to 
h ~ 5.7 x 10~ 5 . Eq. ([6]) with the two values of h are 
shown as dashed lines in figure [5J 

In Ref. [l| a theory was developed to show how the 
abundant e + e~-pairs during the BBN, for temperatures 
below the neutrino decoupling (Tg ~ 0.7 MeV), result 
in modifications in the Salpeter formula. In fact, only a 
very small fraction of the electrons present in the medium 
is needed to neutralize the charge of the protons. The 
major part of the electrons are accompanied by the re- 
spective positrons created via 77 — > e + e~ so that the 
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FIG. 2: (Color online) Baryon density (solid curve) during 
the early universe as a function of the temperature in units 
of billion degrees Kelvin, Tg. The dashed curves are obtained 
from Eq. @ with ft ~ 2.1 x 10~ 5 and ft ~ 5.7 x 10" 5 , 
respectively. 

total charge of the universe is zero. At the decoupling 
temperature the neutron density is only about 1 /6 of the 
total baryon density. Thus, the ion charge density (pro- 
ton density) at this epoch is approximately equal to the 
baryon density. With this assumption the enhancement 
factor in Eq. Q becomes independent of the tempera- 
ture, 

f BBN = exp (4.49 x l(r 8 C^ 2 ) 

- 1 + 4.49 x 10" 8 C^i^2 for T 9 < 1, (7) 

and 

jBBN = exp ( 2 71 x iq-8^ZxZ 2 ) 

~ 1 + 2.71 x KT 8 CZi^2 for T 9 > 2, (8) 

which yield small screening corrections for all known nu- 
clear reactions. 

It is also worthwhile to calculate the Debye radius 
as a function of the temperature. This is shown in 
figure [3] where we plot Eq. Q with the ion density 
equal to the proton density. The accompanying dashed 
lines correspond to the approximation of Eq. with 
h ~ 2.1 x 1CT 5 and h ~ 5.7 x 1(T 5 . This leads to two 
straight lines in a logarithmic plot of 

R D = R^T 9 -\ (9) 

with R$ ~ 6.1 x 1(T 5 cm and ~ 3.7 x 1(T 5 cm, re- 
spectively. In figure [3] we also show the inter- ion distance 
by the lower dashed line. It is clear that the number of 
ions inside the Debye sphere is at least of the order 10 3 , 
which would justify the mean field approximation for the 
ions. In contrast to protons, electrons and positrons are 
mostly relativistic and their chaotic motion will proba- 
bly average out the effect of screening around the ions. 
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FIG. 3: (Color online) Debye radius during the BBN as a 
function of the temperature in units of billion degrees Kelvin 
(solid line). The dotted lines are the approximation given by 
Eq. © with R^ 1 ~6.1x 10" 5 cm and ~ 3.7 x 1CT 5 cm, 
respectively. The inter-ion distance is shown by the isolated 
dashed-line. 



But because the number density of electrons is large, an 
appreciable fraction of them still carry velocities compa- 
rable to those of the ions. The effects of such electrons 
on the modification of the Debye- Hueckel scenario might 
be worthwhile to investigate theoretically. 

Based on the above arguments, the screening by elec- 
trons during the BBN is likely a negligible effect. But 
one needs to verify arguments that sensitive quantities 
such as the Li/H ratio predicted by BBN might be im- 
pacted. This ratio is very small but is one of the major 
problems for the Big Bang predictions. In fact, there are 
discrepancies between the BBN theory and observation 
for the lithium isotopes, 6 Li and 7 Li. These discrepancies 
are substantiated by recent observations of metal-poor 
halo stars [32j and the high precision measurement of 
the baryon-to-photon ratio r\ of the Universe by WMAP 
[3ll |33[ 34] . In view of the relevance of this topic to BBN 
and its predictions, it is worthwhile to check the influence 
of the electron screening on the elemental abundances. 

The BBN is sensitive to certain parameters such as 
the baryon-to-photon ratio, number of neutrino fami- 
lies, and the neutron decay lifetime. We use the values 
7] = 6.19 x 10~ 10 , N v = 3 and r„ = 878.5 s for the 
baryon/photon ratio, number of neutrino families and 
neutron-day lifetime, respectively. We have included all 
reactions of the standard BBN model using the values 
of the cross sections as published in Ref. Q. The re- 
action rates were modified to include screening factors 
calculated with Eq. (|7|). 

In order to test under which circumstances the screen- 
ing by electrons would make an appreciable impact on 
the BBN predictions, we have artificially modified the 
Slapeter formula for f BBN by rescaling it with a fudge 
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FIG. 4: (Color online) Variation (in percent) of the abun- 
dances of several light nuclei as a function of a multiplica- 
tive factor w artificially enhancing Salpeter's formulation of 
screening, w is defined in Eq. (|10p . 



factor w, i.e. 

/' = exp {^ Zl I ^ j J r 1 ■ or inf' = wlnf, (10) 

with w varying between 1 and 10 4 . We quantify the de- 
pendence of BBN on the electron screening by the quan- 
tity 



A = 



Y' -Y 

Y ' 



(11) 



where Y denotes the abundance Y of an element pro- 



duced during the BBN. The primed quantities Y' denote 
the abundances modified by the screening effect. 

In figure [4] we show the value of A (in percent) as a 
function of the screening enhancement factor w in Eq. 
([7]). For w ~ 1 any of the abundances are modified by 
less than 10~ 5 in case that the screening factor is cal- 
culated according to Eq. If the w is increased the 
abundances of 6 Li and 7 Be increase, while those of D, T, 
3 He and 7 Li decrease as w increases. This shows that 
the elemental abundance ratios have a different sensi- 
tivity to the electron screening. These calculations also 
make it evident that the relative abundances in the BBN 
would be somewhat modified only if the electron screen- 
ing would be enormously enhanced (by a factor w larger 
than 10 4 ) compared to the prediction of Salpeter's model. 

In conclusion, using a standard numerical computation 
of the BBN we have shown that electron screening can- 
not be a source of measurable changes in the elemental 
abundance. This is verified by artificially increasing the 
screening obtained by traditional models 11]. We back 
our numerical results with very simple and transparent 
estimates. This is also substantiated by the mean-field 
calculations of screening due to the more abundant free 
e+e~ pairs published in Ref. [l|. They conclude that 
screening due to free pairs might yield a 0.1% change on 
the BBN abundances. But even if mean field models for 
electron screening were not reliable under certain condi- 
tions, which we have discussed thoroughly in the text, it 
is extremely unlikely that electron screening might have 
any influence on the predictions of the standard Big Bang 
nucleosynthesis . 
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